Time-Dependent Transport Through Molecular Junctions 
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We investigate transport properties of molecular junctions under two types of bias — a short time pulse or 
an AC bias — by combining a solution for the Green functions in the time domain with electronic structure 
information coming from ab initio density functional calculations. We find that the short time response depends 
on lead structure, bias voltage, and barrier heights both at the molecule-lead contacts and within molecules. 
Under a low frequency AC bias, the electron flow either tracks or leads the bias signal (capacitive or resistive 
response) depending on whether the junction is perfectly conducting or not. For high frequency, the current lags 
the bias signal due to the kinetic inductance. The transition frequency is an intrinsic property of the junctions. 



The goal of achieving the ultimate miniaturization of elec- 
tronic components is the driving force behind the realization 
of molecular electronic devices. The idea dates back to at least 
1974 I]]], and technology has advanced especially rapidly in 
the last decade |2|.|3|]. While most nanoscale transport studies 
have focused on steady state behavior, recently the high fre- 
quency (GHz or THz) performance of nanotube and graphene 
diodes or transistors has been investigated Willi . The small 
junction areas, low capacitances, and high electron mobilities 
of these molecular devices seem to offer a cutoff frequency in 
the THz range fl3l . From the theoretical point of view, the 
question of how molecules behave under time-dependent per- 
turbation has to be answered since the short time response of 
functional units is essential to the construction of molecular 
electronic devices. 

In recent years, different theoretical approaches have been 
developed for this purpose. The methods and schemes em- 
ployed in these theoretical studies include, for example, time- 
dependent density functional theory (TDDFT) combined with 
a Green function technique for open model systems lfl3ll . 
a nonequilibrium Green function (NEGF) method treating 
the time-dependent bias as a perturbation to the steady-state 
Hamiltonian 11411 . real-time TDDFT propagation for closed 
systems fl5l llfJl . a quantum master equation scheme based 
on TDDFT for model systems f 1711 . a combination of TDDFT 
and NEGF with the wide-band limit approximation for open 
systems lfl8l[l9"tl . real-time propagation of the Kadanoff-Baym 
equations for open and interacting model systems [20-22], 
and a self-consistent NEGF formalism for the electron trans- 
port through nanotubes under a time-dependent gate potential 
112311 . Despite the large theoretical effort made by different re- 
search groups, the computational studies for real open systems 
with atomic details described by ab initio electronic structure 
calculations lfl8l [l9l l23ll are still lacking because of the com- 
putational difficulty. 

Within the Keldysh nonequilibrium Green function descrip- 
tion l24ll . Jauho and coworkers have derived formulations of 
transport in the mesoscopic regime under influence of external 
time-dependent perturbations. Based on this work, Zhu and 
coworkers [ 14] established a computationally efficient method 
without the need of the wide-band limit (WBL) approxima- 
tion by using the zero bias equilibrium Green function as the 
initial state of a tight-binding model system and carrying out 



the analysis in the time domain. The finite correlation time 
in open systems reduces in large part the computational effort 
and makes this method a practical approach for ab initio study. 

In the present work, we follow this idea and extend the 
tight-binding level theory to a density functional theory (DFT) 
description of the electronic structure, i.e., the initial equilib- 
rium states are obtained by a DFT-Green function formula- 
tion dHHH. 

We investigate a relatively simple atomic chain 
system so that we can conveniently adjust various parame- 
ters in order to explore the general behavior of molecular 
junctions under a time-dependent bias, including a short time 
pulse and an AC bias. In this way, we discuss how the time- 
dependent transport is affected by the nature of the leads, the 
lead-molecule coupling, the barrier in the molecule, and the 
amplitude and frequency of the bias signal. 

Following Ref. IU4l 12411 . the time dependent current is 

/«(*) = y^ eTr / dt 1 [G r (t,t 1 )X<(t 1 ,t)+G<(t,t 1 )-Z a a (t 1 ,t)}, 

(1) 

where E < ' a (ti,t) is the lesser/advanced self-energy. When 
no time-dependence is present, the steady state Green func- 
tion G is solved by using the NEGF technique in energy space 
ll29ll . Under a time-dependent voltage, the single-particle en- 
ergies become time-dependent in the leads, which causes ac- 
cumulation and depletion of charges to form a dipole across 
the device region. Based on the known steady-state Green 
function G, we solve for the retarded lesser Green functions 
by the Dyson and Keldysh equations, respectively lfl4l l24ll . 

We first investigate the effect of lead structure on the tran- 
sient response. The system studied is a diatomic hydrogen 
molecule (H2) with a 1 A bond length sandwiched between 
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FIG. 1: Illustration of a diatomic hydrogen molecule (yellow) with 
1 A bond length sandwiched between hydrogen chain leads. We vary 
the interatomic distance in the hydrogen chain but keep the distance 
between the H2 and the leads fixed, d = 1.5 A. The H2 molecule 
together with four hydrogen atoms on each side form the extended 
molecule (in the blue box). Bias is applied to the left lead. 
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FIG. 2: Current as a function of time (both in atomic units) for the H2 
molecule of Fig.[TJand three different leads: the interatomic distance 
in the leads is (a) 1.5 A, (b) 2.0 A and (c) 2.5 A. A square shaped 
voltage pulse is applied to the left lead starting at t — and ending at 
425. Solid, dotted, and dashed lines correspond to applied voltages 
of 1, 3, and 5 V, respectively. 

two semi-infinite hydrogen chain leads as shown in Fig. Q] 
The distance from the leads to the H2 molecule is fixed to be 
1.5 A. We change the H-H distance in the leads, thus vary- 
ing the interatomic coupling strength. The narrow band width 
produced by a weak interatomic coupling constrains the elec- 
trons, producing less metallic leads. 

Single-zeta basis sets (SZ) and optimized Troullier-Martins 
pseudopotentials |f3~b]l are used for all calculations lHH]. For 
hydrogen, if the applied voltage is low so that p states are not 
excited, a SZ basis set is a good approximation. The PBE ver- 
sion of the generalized gradient approximation (GGA) func- 
tional [32] is adopted for exchange-correlation. For conve- 
nience, we adopt atomic units for electric current and time (e 
= h = m e =1). An accurate description of atomic and chemical 
details in the contact region obtained from DFT calculations 
enables a meaningful study of realistic molecular devices. 

The I(t) characteristics for leads with different H-H dis- 
tances are shown in Fig. [2] A square shaped voltage pulse is 
applied to the left lead starting at t = (a.u.) and ending at 
t = 425. The time mesh varied between 200 and 600 points 
for all results shown here; generally, more points were needed 
for more strongly coupled leads. 

The current initially increases for a short time and then os- 
cillates because of interference for several to tens of a.u. de- 
pending on the type of leads. We call the time needed for the 
current to increase from zero to its maximum the finite cur- 
rent response time, tr. This delay in response is related to 
the inertia (effective mass) of the carrier and provides a mech- 
anism for inductance (kinetic inductance) different from the 
magnetic one. As tr is very short (of order 1 fs), the kinetic 
inductance is only observable at high frequencies. 

The second feature to note in Fig.|2]is that the initial rise in 
current is steeper when the interatomic distance in the leads 
is smaller (smaller tr) for the same applied voltage. For in- 
stance, 1 V bias, tr is about 4, 10, and 25 for interatomic dis- 



-6 -1-2 2 -1 IS 

FIG. 3: Transmission as a function of energy for the H2 molecule of 
Fig. [Hand the same cases as in Fig.[2] Solid, dotted, and dashed lines 
correspond to applied voltages of 1, 3, and 5 V, respectively. 

tances 1.5,2.0, and 2.5 A, respectively. This behavior is easily 
understood by considering the effective mass of the electrons 
in these different leads: as the interatomic distance increases, 
the band width, of course, decreases, and so the electrons have 
a larger effective mass. 

Fig.|2]also shows the relations between bias voltage and the 
response time. For the leads with wider bands (band width 
wider than the bias window), tr is almost independent of the 
bias (panels a and b), indicating that the kinetic inductance is 
basically a constant. On the other hand, for the narrow band, 
less metallic lead (panel c), tr decreases rapidly for larger 
bias, showing that the kinetic inductance in this case depends 
strongly on the magnitude of the bias. One contributing fac- 
tor is that the electrons are accelerating in the electric field; 
another is the mismatch between the narrow band in energy 
on left and that on the right when the bias window is larger 
than the band width (see Fig. [3] for further discussion). The 
light might occur, for instance, if a narrow band metallic ox- 
ide or silicide is used as the lead material. Over the duration of 
the transient response, the bias causes alternately accumulat- 
ing and depleting regions of charge. Subsequently, the current 
finally tends to a steady state. 

In our system, a higher bias does not necessarily lead to a 
larger steady state current — negative differential conductance 
can occur. In fact, the possibility of highly nonlinear I-V 
curves is one of the appealing features of molecular electron- 
ics. Fig. [3] shows the transmission functions, T(E), for the 
same cases as Fig. [2] For reference, the widths of the s-band in 
infinite hydrogen chains with interatomic distances 1.5 A, 2 A 
and 2.5 A are 12.8 eV, 5.7 eV, and 2.4 eV, respectively. Focus- 
ing on the case with the narrowest band in the leads (2.5 A 
interatomic separation), we see that the net current becomes 
nearly zero for a bias of 3 V (Fig. [3};). This is because the 
s-bands for the left and right leads have no overlap when the 
bias is greater than the band width, 2.4 eV. Incident electrons 
encounter a hard wall, causing the current to oscillate over a 
much longer time than for wide bands. More generally, we 
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FIG. 4: Schematic of two extended-molecules: (a) H2 and (b) H10 
shown with four atoms of each lead, dl and d2 are the lead-molecule 
distances. 

see that the width of the transmission window equals the dif- 
ference between the s-band width and the applied bias. 

The discussion in the last paragraph suggests the question: 
Does a less transparent junction generally lead to a longer 
time to reach a steady state? To answer this, we vary the 
barrier height associated with the molecule while keeping all 
other parameters constant. The barrier height at the molecule- 
lead contacts can be altered easily by changing the distance 
between the molecule and the leads. The barrier presented 
by the molecule itself can be changed by varying the num- 
ber of atoms. Thus we shall compare results for H2 and H10 
molecules (Fig.|4]l. 

For isolated H2 and H10 with an interatomic spacing of 
lA, the gaps between the highest occupied molecular or- 
bital (HOMO) and the lowest unoccupied molecular orbital 
(LUMO) are 10.7 eV and 3.2 eV, respectively (DFT GGA cal- 
culation). Fig. 2] shows only the extended molecules; note 
that the bias voltage is applied in the leads, i.e., outside of the 
extended-molecule region. 

Fig. [5] shows the I(t) curves of both molecules for three 
values of the lead-molecule separation. The bias voltage ap- 
plied at t — is not turned off during the simulation. When 
the molecule-lead coupling is strong (dl = d2 = 1.5 A, panel 
a), both the H2 and H10 junctions reach steady state quickly, 
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but note that higher bias causes a longer transient regime. The 
magnitude of the steady state current is nearly the same for 
the H2 and H10 molecules. 

When the molecule-lead distance is stretched by 0.2 A to 
1.7 A (panel b), the low bias behavior of the molecules re- 
mains virtually the same, and further there is only modest 
quantitative change from the strongly coupled case. However, 
for a 3 V bias, the I(t) curves differ substantially, both from 
each other and from the previous case. In the short molecule, 
the duration of the overshoot and oscillating regime is greatly 
extended, while in the long molecule, there is a surprisingly 
small decrease in the steady state current (2.84xl0 -2 in panel 
a to 2.25 x 10 -2 in b). H2 experiences a more substantial 
change in current (from 2.70 x 10~ 2 to 1.65 x 10~ 2 ) in line 
with that expected from the fractional change in current at low 
bias. Thus, the higher molecule-lead barriers are more clearly 
manifest in the short molecule and hidden in the long one. 

These differences are amplified further in the case of largest 
molecule-lead distance (2 A, panel c). The low bias traces are 
quite similar to each other. At high bias, I(t) for H2 oscil- 
lates for a long time, while the steady state current for H10 is 
surprisingly high. 

We believe these difference between H2 and H10 are caused 
by the larger HOMO-LUMO gap in the short molecule. Fig. [6] 
shows the projected density of states and transmission func- 
tions of both molecules in the weak coupling case (molecule- 
lead separation of 2 A). The T(E) within the bias window for 
a 1 V bias (i.e. —0.5 V to 0.5 V) are almost identical for H2 
and H 10 . In contrast, for a 3 V bias, the T(E) within the bias 
window, now —1.5 V to 1.5 V, are totally different: for H 10 , 
the tails of two resonant peaks extend into the bias window. 
Therefore the long molecule is much more transparent, lead- 
ing to a larger current flow and less oscillatory behavior com- 
pared to the short molecule. Thus we see that both molecule- 
lead barriers and internal barriers within the molecules can 
cause significant changes in the I(t) characteristics. 

After studying the transient response to a square shaped 
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FIG. 5: Current as a function of time for both H2 and H10 molecules 
with different coupling strengths to the leads: the molecule-lead dis- 
tance is (a) 1 . 5 A, (b) 1 . 7 A, and (c) 2 A. Results for a bias of both 1 V 
(H2 solid, and H10 dashed) and 3 V (H2 dotted, and H10 dot-dashed) 
are shown. The bias is applied at t — 0; the interatomic distance in 
the leads is 1.5 A. 



FIG. 6: Projected density of states and transmission functions for 
both the long and short molecules connected to hydrogen chain leads 
(1.5 A interatomic distance). The molecule-lead distance is fixed at 
2 A. At energies more than 1 eV away from Ef, the properties of H2 
molecular junctions (IV solid, 3 V dotted) differ greatly from those 
of H10 (1 V dashed, 3 V dot-dashed). 
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FIG. 7: Current as a function of time for H2 molecules subjected to 
AC signals of varying frequencies, smaller to higher from panel a to 
e. A well-coupled case (solid, dl — 1.5 A) and weakly-coupled case 
(dashed, dl = 2 A) are compared to the applied voltage (dotted, 1 V 
amplitude). The interatomic distance in the leads is 1.5 A. 

pulse, we now look at how a molecular junction acts when 
a sinusoidal voltage with period T is applied. The frequency 
should be slower than the plasma frequency uip of the leads 
so that the electric field is effectively screened and the voltage 
drops across the device region. For most metals, the plasma 
frequency is in the ultraviolet regime, ranging from 10 15 Hz 
to 10 17 Hz, so typically this criterion is satisfied. For the ID 
hydrogen chains in our calculations, is of order 10 17 Hz, 
assuming the dielectric constant and permeability are that of 
vacuum and the length of the device is around 1 nm lf33fl . 

Fig.|7]shows the AC response of H2 junctions (schematic in 
Fig- [D with two molecule-lead distances, 1.5 and 2 A, repre- 
senting the well-coupled and weakly-coupled regime, respec- 
tively. First, consider the low frequency cases, T/4 tr 
(panels a, b, and c). The current response of the well-coupled 
system tracks the sinusoidal signal, while for the weakly cou- 
pled junction, the current leads the bias voltage. This behavior 
can be understood by considering the equivalent electric cir- 
cuit of the lead-molecule-lead system. At low frequency, the 
inductance is not important, so we ignore it temporarily; then 
we can view the junctions as a resistor (R) and capacitor (C) 
in parallel formed by the two contact interfaces. For the well- 
coupled system (dl = 1.5 A), there is actually no interface, 
and thus the capacitor disappears, and the whole system is ba- 
sically resistive. When the coupling is weaker, the capacitor 
is formed, making the system capacitive in nature. 

As the frequency increases, the kinetic inductance grad- 
ually appears; consequently, the phase shift of the weakly- 
coupled system becomes progressively smaller in panels a, b, 
and c of Fig. [7] As the frequency increases further, T/4 be- 
comes shorter, approaching the finite current response time 
tr. Then the effect of the delay in current becomes signifi- 
cant and at a certain frequency the characteristics of the junc- 
tion changes from capacitive to inductive, as shown in panel 
d. Our results are qualitatively consistent with a recent cal- 
culation for an Al-nanotube-Al junction where it was found 



mi im> 1 an m> :\id 4m) 4,-ri) 



-o.u: 



iJAV Jtv'AV ztVvr 



100 150 200 250 300 .550 40O 150 500 

L (il.U.) 



FIG. 8: Current as a function of time (solid) for H2 molecules sub- 
jected to AC signals (dashed) of varying frequencies, smaller to 
higher from panel a to f. The molecule-lead distance dl equals 1 .5 A, 
and the interatomic distance in the leads is 2.5 A. 

that the system is inductive for a high bias frequency (T ~ 
response time) lfl9ll . 

To support our explanation, we consider different values of 
tr. Since we showed that tr in ID molecular junctions is 
significantly affected by the nature of the leads, we consider a 
different lead structure but keep the intramolecular and lead- 
molecule structure the same. Thus the lead-molecule separa- 
tion (1.5 A) becomes different from the interatomic distance 
(2.5 A) in the leads, and therefore a capacitor is formed. The 
I{t) curves in Fig.[8]show that the current now leads the bias 
even for very small frequency (panels a and b). The current re- 
sponse changes, as expected, from capacitive to inductive as 
the bias frequency increases; however, compared to the sys- 
tem with 1 .5 A interatomic distance, the transition frequency 
decreases from 2.8 x 1CT 2 (Fig. 01) to 3.6 x 1(T 3 (Fig.©). 
This reduction is consistent with the change in tr\ As the in- 
teratomic distance in the leads becomes larger, tr increases 
from ~ 4 to ~ 25 (Fig. [7] a and c). As a caution, we men- 
tion that since the intrinsic magnetic inductance has not been 
taken into account here, the total inductance of the ID molec- 
ular junctions may be underestimated; however, the magnetic 
inductance may be much smaller than the kinetic one for high 
frequencies (THz) fi - 

It is interesting to note that the magnitude of the current 
driven by an AC bias is not necessarily consistent with that 
driven by a DC bias. In Fig. |5J the current through H2 with 
dl = 2.0A is smaller than that with dl = 1.5A. This is because 
the contact barrier in the former is higher. In contrast, under a 
high-frequency AC bias, the magnitudes of the current in these 
two cases become very close, as shown in Fig. [7] Although 
both bias frequencies are much slower than the plasma fre- 
quency uip, the time needed for the system to reach its steady 
state is much longer, by orders of magnitude, than 1 /top due 
to the very weak screening in the 1-D structure. Only when 
T/4 of the AC bias is much longer than this time does the 
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FIG. 9: Current as a function of time for H2 molecules subjected to 
AC signals of different magnitude: 1 V, 3 V, and 5 V for panels a, b, 
and c, respectively. I(t) for two lead-molecule spacings are shown, 
dl = 1.5 A (solid) and 2 A (dashed), and compared to the applied 
voltage (dotted). The interatomic distance in the leads is 1.5 A. 

AC current reflect the steady state current; otherwise, it will 
largely be determined by the transient states. Because the tran- 
sient states are determined by the extended molecule instead 
of the molecule and the contact alone, the magnitude of the 



AC current can be very different from that of the DC current 
in the same system. 

Finally, we point out that when the AC frequency is rela- 
tively small, irregular behavior can be significant, as we saw 
for pulsed signals. We show in Fig.|9]that as the voltage mag- 
nitude increases (from 1 V to 5 V), the current response may 
not follow a sinusoidal wave. The irregular features are larger 
when the junction is less transparent. This is in line with our 
previous conclusion (Fig. [3J;) and applies to the situation of 
relative low frequency AC response. 

In conclusion, the transient response of a molecular junc- 
tion depends on the lead structure, bias voltage, and barrier 
height seen by the transported electrons. A higher electron 
density or a smaller effective mass leads to a faster response 
characterized by a smaller current response time tr. A high 
barrier height yields long oscillatory behavior in current, seen 
in both the pulsed and AC situations. The current follows the 
AC signal only when a junction is perfectly conducting and 
the signal frequency is slow (T/4 tr); otherwise, a lead- 
molecule-lead junction should be viewed as a complex circuit 
consisting resistors, capacitors, and inductors. Currents can 
lead or lag the AC signal, and the transition frequency be- 
tween the two regimes is an intrinsic property of the junction. 
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